Programmable quantum state discriminator by Nuclear Magnetic Resonance 
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In this paper a programmable quantum state discriminator is implemented by using nuclear 
magnetic resonance. We use a two qubit spin-1/2 system, one for the data qubit and one for 
the ancilla (programme) qubit. This device does the unambiguous (error free) discrimination of 
pair of states of the data qubit that are symmetrically located about a fixed state. The device is 
used to discriminate both, linearly polarized states and elliptically polarized states. The maximum 
probability of the successful discrimination is achieved by suitably preparing the ancilla qubit. It 
is also shown that, the probability of discrimination depends on angle of unitary operator of the 
protocol and ellipticity of the data qubit state. 

I. INTRODUCTION 

Researchers have studied the possibility of performing computations using quantum systems and conjectured that 



^ a machine based on quantum mechanical principles might be able to solve certain types of problems more efficiently 
FTithan can be done on conventional compntersfl fl Q. Later Llovd proposed that snch a q »a,t»,n compute, might 
be built from an array of coupled two state quantum systems U|. Its theoretical possibility has generated a lot of 



■ be built from an array or coupled two state quantum sys ti 
$H 'enthusiasm for its experimental realization!?! la M. Isl Ift llo| ■ 

a ■ 



J . In parallel with quantum computation, the related field 
of quantum information theory is developed, which forms the quantum analogue of classical information theory 
Several techniques are b eing exp loited for quan tum computing and quantum information processing, including nuclear 
magnetic resonance [id, La, L^, La La, Lll La Ld| 

Recently quantum state discrimination has been studied extensively in the context of quantum communication and 
quantum cryptography [20I . I21I l^ . I23I I24I l^l . I^rl [^J ■ Quantum state discrimination is the problem of determining 
the quntum state, given the constraint that it belongs to the previously specified set of non-orthogonal states. One of 
the characteristic features of quantum mechanics is that, it is impossible to devise a measurement that can distinguish 
non-orthogonal states perfectly |lfj |. However one can distinguish them with a finite probability by appropriate 
measurement strategy. There are two different optimal strategies of discrimination: (i) Probabilistic discrimination 
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(conclusive result, but error may appear) and (ii) Unambiguous discrimination(inconclusive result may appear, but 
no error). The unambiguous discrimination of two pure states was investigated by Ivanovic |2(|. and the optimal 
procedure was given by Peres [2^ . Cheffels and Barnett have generalized Peres's solution to an arbitrary number of 
equally probable states which are related by a symmetry transformation 28J . The first experiment to discriminate two 
non-orthogonally polarized single photons states of light was done by Huttner et.alj^J. 

Quantum measurement is the final step of any quantum computation. In many situations the choice of an optimal 
measurement depends on the task to be performed. A quantum multimeter is a quantum measurement device which 
can perform a specific class of generalized measurements in such a way that each member of this class is selected by 
a particular quantum state of a programme register |3Cl I'ilL 1.321 l8.il 1.3- 1| . The parameters determining the character 
of quantum measurement can be encoded in a quantum state of a programme register |2 EmI 13 • Dusek et. al [3^] 
have shown that pair of non-orthogonal states of a qubit can be discriminated and the measurement to be done for 
this discrimination is decided by the state of theprogramme qubit (ancilla qubit). One can discriminate several pairs 
of states of a qubit by using the same protocol [3JJ. Such a quantum device is known as quantum multimeter for the 
discrimination of pair of qubit states. Recently Dusek et. al 38] have also demonstrated experimentally the possibility 
to control the discrimination process by the quantum state of ancilla qubit, in linear optics by performing the partial 
measurement in the Bell basis. C rypt ographic applications of quantum state discrimination have been extensively 
studied in the literature [3, QQQ 

Nuclear magnetic resonance (NMR) has played a leading role for practical demonstration of quantum algorithms 
and gates [12I lid I141 1 1.4 llfii livllisi Il4 ] . . The unitary operators needed for implementation of these quantum circuits 
have mostly been realized using spin selective as well as transition selective radio fr equency pu lses and coupling 
evolution, utilizing spin-spin (J) or dipolar couplings among the spins |l2ul3lll4Lll5lll6lll7Ul8Ul9l| . In this paper we 
demonstrate the implementation of quantum state discriminator which discriminates the pair of non orthogonal states 
as well as orthogonal states which are symmetric about a particular state, conditioned on the state of the ancilla qubit. 
We use spin selective pulses and evolution under J-coupling for the implementation. Projective measurement required 
for the discrimination is simulated by a method given by CollinsQ- 

Our experimental results are in agreement with 
the theoretical results 30|. To the best of our knowledge this is the first experimental demonstration of programmable 
quantum state discriminator by NMR. 

In section(II), we discuss the theory of discrimination of both elliptically and linearly polarized states. Experimental 
details and results of different experiments are given in section(III). Results are concluded in section(IV). In the 
Appendix, unitary operators of ideal pulses are derived. 
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II. THEORY 

The following protocol discriminates pair of elliptically polarized states of the data qubit unambiguously (error 
free). Let the two states and \ip2) of the data qubit be (fig. la), 

\ipi) = (xcosQi + iysinB i)\0e>) + (xsin9i — iyco89i)\ljj), 

1^2) = {xcosOi + iysin6i)\0u) ~ (xsinBi — iycos9i)\ljj) . (1) 

Ellipticity (e) of the data qubit states is defined as, tan(e)=y/x. Here y=x corresponds to circularly polarized states 
and y=0 (e= 0) corresponds to linearly polarized states (fig. lb). The protocol uses one ancilla (programme) qubit 
for the discrimination. A quantum circuit for the discrimination is shown in Fig. 2. In this circuit the first qubit is 
the data qubit (I^d)) and the second qubit is the ancilla qubit (\ij)A))- The data qubit can be either or \ip2}- The 
aim of the protocol (fig. 2) is to determine whether the data qubit is \tpi) or \ip2), knowing the angle between 
and |^>2)- It is shown that, a pair of data qubit states |^i) and \ip2) can be discriminated by suitably preparing the 
ancilla qubit. One can switch the apparatus to work with several different pairs of data qubit states. In this paper 
we experimentally demonstrate the discrimination of both elliptically and linearly polarized states, and compare the 
results with simulations. In the following, we first describe the discrimination of elliptically polarized states and later 
the linearly polarized states. 

Elliptically polarized states l^i) and ^2) (eqn. 1) can be re-written as, 

|^i)=o 1 |0)+6i|l), 

\ih) = ai\0) - h\l), (2) 

where a\ = xcosQ\ + iysinOi and b\ = xsin6\ — iycos9\ are complex numbers, where by definition x 2 + y = 1. 
ai and 61 can also be written in polar form as, 

a\ = e 1 ^ 1 cos{r)) and b\ — e 1 ^ 2 sin(i]), where, (3) 



tan{n) = p-r, tan{4> x ) = ( ^"f 1 ), and tan(0 2 ) = ( VC ° S Q 1 )- (4) 
ai xcosui xsmui 



Then and \1jj2} can be written as general states on the Bloch sphere 



|^) = cos(r))\0) + e*«n(»7)|l), 

|^ 2 )=co«(» J )|0)-e i *«n(»/)|l), (5) 



where 4> = 4> 2 — 4>i with the overall phase (e 1 ^ 1 ) being neglected. 
Let the ancilla qubit (programme qubit) be, 

\iI>a) = a 2 \0) + b 2 \l), 



(6) 



where a 2 = xcos6 2 + iysin9 2 and b 2 = xsin9 2 — iycos6 2 . To discriminate and \ip 2 ), the condition on a 2 and 62 is 
derived as follows. 

The total input state is \iPda) = \^d) ® \iPa), where the data qubit \ipn) is cither or \ijj 2 ) (eqn. 2). 



DA 



{ ai a 2 \0 D A ) + ai b 2 \0 D l A )) ± (b ia2 \l D A ) + 6i6 2 |1bU)). 



(7) 



The sign of the second term in Eqn. (7) determines whether the data qubit is l^i) or \ip 2 ). The protocol for the 
discrimination requires a unitary transformation on \iPda) given by |3(| . 



U 



(8) 



/ cos(a) —sin(a) 0\ 
sin(a) cos(a) 
10 

V 01/ 

where a is a fixed parameter which does not depend on the data and programme qubits states. The unitary 

operator given in Eqn. 8, is a rotation in the subspace spanned by lO^O^) and |0d1a), which is achieved here by two 
controlled- not gates and four single qubit gates (fig. 2), where the Single qubit gates are given by, 



cos(a/2) sin(a/2) 
—sin(a/2) cos(a/2) 

And controlled- not (CNOT) gate is given by, 









*)•«-{ 







,u2 = 



cos(a/2) —sin(a/2) 
sin(a/2) cos(a/2) 



(9) 



CNOT 



/l 0\ 
10 



1 

Vo 1 0/ 

After the application of the unitary transformation U(eqn. 8), the final state is, 



(10) 



U\^da) = {a\a 2 cosa ~ aib 2 sina)\Q D QA) + {a\a 2 sina + aib 2 cosa)\QD^A)^ I&i^IIdOa) + bib 2 \\D^A)\- (11) 
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For successful discrimination of data qubit states \ip\) and 1^2), the condition on the coefficients of Eqn. (11) is, 

{a\a 2 cosa — aib 2 sina) = 6102, (12) 

This yields, 

U\ipDA) — [6iO2|0_d0 j 4) + (aia 2 sina + aib 2 cosa)\0DlA)} 

± [6io 2 |1d0a) +6i62|1dU>]. (13) 

Equation (13) can be re-written as, 

U\iPda) = V2bia 2 \±) + (ciia 2 sina + aib 2 cosa)\0D^A) ± &i&2|1.d1a), (14) 

where |±) = 1/V2(|0 D A ) ± \1 D A )). 

If the final state (eqn. 14) contains state |+) then the initial state of the data qubit is \ipi) if the final state contains 

state |— ) then the initial state of the data qubit is \ip 2 ). Square of the coefficient (^/2bia 2 ) of the state |±) is called 

the probability (P) of discrimination [3^. 

Equation (12) gives the condition on the ancilla qubit state. For example, when a = 90°, (b 2 /a 2 ) = — (bi/a-i) , and 

from Eqn.s 2 and 6 it is seen that \ipA)=\il>2}. However for other values of a, \4>a) differs from |?/>i) and {ifa}. Hence 

a = 90° is a special case of Eqn. (12). 

In NMR, the measurement is performed on an ensemble and the results are contained in the expectation values 

((<J X ) and (<T y )) of Pauli spin matrices, which in the frequency space yield intensities of various transitions. From Eqn. 

(14) it is noted that the two transitions of the data qubit have different intensities. The lO^O^) *-* |1d0a) transition 

has the intensity ±6^02, and the |0z3 lyi) <-> |1d1a) transition has the intensity ±(aia 2 sin(a) + aib 2 cos(a))bib 2 . To 

find whether the final state (eqn. 14) contains |+) or |— ), one has to do the projective measurement on the state 
±). To simulate projective measurement, we use a method given by Collins for an expectation value quantum search 
43l |44j. In our experiments, the goal of the projective measurement is to collapse the state of the ancilla qubit 

(given by equation 14) to |0a) so that the data qubit gives only one peak which is the coherence of the superposition 

state |±) (or IOdO^) <-> |1d0a) transition). One can collapse the state of the ancilla qubit to \0a) by adding two 

experiments (detection on the data qubit), one without controlled-er 2 gate and one with controlled-cr 2 gate (fig. 2). 

Here the controlled-cr z gate is given by the unitary transformation, 



/l \ 

10 

10 

Vo -1/ 



(15) 



G 



Controlled-<7 z gate inverts the sign of the statell^d^)- When <j c z is applied after the unitary transformation U (eqn. 
8), the final state (eqn. 14) becomes, 

(tIU\iPda) = V26ia 2 |±) + (aia 2 sina + aib 2 cosa)\QD^A) T bib 2 \lD^A)- (16) 

In the first experiment (eqn. 14), intensities of the data qubit transitions corresponding to |0 zjO^) <-> |1£>0a) and 
|0_d1a) |1_d1a) are, ±b\a^ and ±(aia,2sin(a) + ai&2Cos(a))&i&2 respectively, whereas in the second experiment 
(eqn. 16) intensities are ztb^a^ and z f(aia2sin(a) + ai&2Cos(a))&i&2 respectively. Thus when the two experiments are 
added, intensity of IOdIa) I Id 1a) transition goes to zero and that of |0d0a) |1z>0a) to ±2b\a\- Hence by the 
above procedure the ancilla qubit state is collapsed to \Qa), and the phase of the observed transition yields the result 
of the measurement. If the phase is positive then the data qubit is and if it is negative then the data qubit is 
\ip2}- The resultant intensity (26f a^) gives the probability of successful discrimination. 

For the case of linearly polarized states (eqn. 2, y=0; a\ = cos9i, b\ = sin9i), the data qubit states \ipi) and 
^2) are schematically shown in Fig. (lb), and rj and <f> of Eqn. (5) are respectively given by 9\ and zero. In this 
case ancilla qubit state (fig. lb) is given by Eqn. (6), with a 2 — cos02 and 62 = sin6 2 . Rest of the procedure to 
discriminate |"0i) and \ip 2 ) remains the same and the probability of discrimination is given by 2b\a^. 



III. EXPERIMENT 



In NMR spin- 1/2 nuclei having sufficiently different Larmor frequencies and weakly coupled to each other by indirect 
exchange (J) couplings are used as qubits. The Hamiltonian of the two weakly coupled spin-1/2 nuclei is of the form, 

H = U>il zl + U) 2 I Z 2 + 2irJi2l z iI Z 2- (17) 

We have used a Carbon-13 labeled 13 CHCl3 as a two qubit system, where the proton ( 1 ff) and the labeled carbon( 13 C) 
act as two individual qubits. J-coupling between 13 C and 1 H is 209 Hz. The measured longitudinal (Ti) and transverse 
(T 2 ) relaxation times of X B and 13 C are: 1 H (T x =4.8s and T 2 =3.3s), and 13 C (2i=17.2s and T 2 =0.35s). To implement 
the circuit of Fig. (2), the data ( 1 H) and ancilla ( 13 C) qubits have to be first prepared in a pure state. However in 
NMR pure states are difficult to prepare, instead we prepare pseudo pur e states which mimics the pure states. Several 
methods are known for the preparation of pseudopure states 46J, 1471 1481 1491 l50l l5ll |52| . Here we use spatial averaging 
method |53j to prepare pseudopure state using the pulse sequence given in Fig. (3). This pulse sequence |53| is 
specific to labeled 13 C- 1 i? system and different from homo nuclear case. The details of the preparation of pseudopure 
state are given in figure captions. Spectra of equilibrium state and pseudopure state are shown in Fig. (4). After 
preparation of pseudopure state, the quantum circuit of Fig. (2) is implemented by the pulse sequence given in Fig. 
(5). The pulse sequence in Fig. (5) consists of three parts, 
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(i) Preparation of initial state {\4>da))'- After preparation of pseudopure statedOcO^)), the data qubit Q~H) is 
prepared in elliptically polarized state \iPd) (eqn. 5) by applying a 2r\ pulse of appropriate phase on the data qubit 
state |0d)- To prepare the data qubit in state \4>i) or l^) , the phase of the 2-q pulse is (tt/2 + (j>) or —(tt/2 + (f>) 
respectively (Appendix). The ancilla qubit ( 13 C) is prepared in state \i/>a) (eqn. 6) by using Eqn. (12). For example 
for a — 90°, since \ipA) = \ip2}, the ancilla qubit \ipA) is prepared by another (2^)_( x / 2 +0) pulse. For arbitrary a, the 
ancilla qubit is prepared using Eqn. (12) with appropriate pulse angle and phase. In case of linearly polarized state 
(eqn. 2, y = 0), r/=9i and <j) = 0°, the data qubit can be prepared in states and \ip2) respectively by applying 
(29i)y and (29\)- y pulses (Appendix) on \0d}- The ancilla qubit (eqn. 6) is then prepared by applying (26*2)3/ pulse 
on \0a), where 26*2 is calculated according to Eqn. (12). From Eqn. (12), 26*2 can take positive as well as negative 
values. For example When a= 60°, and 20 x = 20°, 40°, 60°, 80°, 90°, 100°, 120°, 140°, 160°, 26> 2 takes the values, 
41°, 17.8°, -10.2°, -42.8°, -60°, -77.2°, -109.8°, -137.8°, -161° respectively. Here one should note that (-29 2 ) y 
pulse is identical to (262) -y pulsc. 

(ii) Applying unitary operator U : The unitary operator U (fig. 2) is prepared by using two CNOT gates, two NOT 
gates and two other single qubit gates u\ and 112 (fig- 5). The NOT gates on the data qubit ^H) are implemented by 
(tt) x pulse and Ui, u 2 on ancilla qubit by (ot)- y and (a) y pulses respectively on 13 C The CNOT gate is implemented 
by using the pulse sequence {t7 /2)\-{-k /2)^-{1/2J)-{-k /2) x -{tt /2) 2 _ z |46j, where the superscript 1 stands for proton and 
2 stands for carbon. The (ir/2)l is obtained by the composite pulse (7r/2)^-(7r/2)^-(7r/2)L 1/ as shown in Fig. 5. The 
(tt/2) 2 _ z pulse can be obtained by a another composite pulse {it /2) 2 _ x -{ir /2)^-{-k /2) 2 , so that the first (n/2) 2 _ x pulse 
of the composite pulse cancels the last (ir/2) 2 , pulse of the CNOT gate yielding the last two pulses in the CNOT 
sequence as (tt / 2) 2 i y - (it/ 2) x . All the pulses in the pulse sequence are applied at resonance, so the chemical shifts 
are refocuscd throughout the pulse sequence. Hence during the time period (1/2J), system evolves only under the 
J-coupling Hamiltonian Hj — 2irJI z iI Z 2 yielding the unitary operator, e _ITr/il/z2 . 

(iii) Controlled-<7 z gate (a z ) is implemented by (7r/2) z ' 2 pulse followed by an evolution for the time 1/2J [l^ . 
(7r/2) z ' 2 pulses are realized by composite rotation on both qubits as shown in Fig. 5. 

(a) Linearly polarized states: 

We have studied the linearly polarized case by varying both the parameters a (rotation angle of U, eqn. 8) and 
29\ (angle between \i\>\) and l^a)) n S- lb). The pulse sequence given in Fig. (5) is implemented, with the initial 
state prepared as described above (in section III (i) ) . Experiment is performed to discriminate several pairs of linearly 
polarized states for a=30°, 45°, 60°, and 90°. For each value of a, the experiment is carried out for 26*1=20°, 40°, 
60°, 80°, 90°, 100°, 120°, 140°, 160°. As mentioned in theory section (II), the experiment is performed twice, one 
with <j c z and other without <r c , and the results are added so that the resultant intensity of the data qubit transition 
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gives the probability of discrimination (Y=2b\a\). Figure (6) contains typical spectra for 26q = 90° and a = 90°, 60°, 
45°, 30°, where the data qubit is prepared respectively in states (fig. 6a-d) and \tp2) (fig- 6e-h). As shown in 
Fig. (6), the positive intensities of the resultant peaks indicate that the initial state of data qubit is l^i) and the 
negative intensities of the resultant peaks indicate that the initial state of data qubit is IV^)- The intensity of the peak 
yields the probability (P=26fa^). In Fig. (6) one can observe that the intensity of the resultant peak (probability 
of discrimination) changes with a. For different values of a, the probability of discrimination P (experimental and 
simulation results) as a function of 29i is given in Fig. (7). From Fig. 7 one can find the optimum angle 29i for 
maximum probability of discrimination for a given value of a. Figure (8), on the other hand, shows the variation of 
Probability of discrimination(P) as a function of a, for different 2B\. From Fig. (8), one can find the value of a to 
get the maximum probability of discrimination for a given angle (2#i) between and \ip2}- In both figures 7 and 
8, the experimental points agree well with the simulations, confirming successful discrimination of linearly polarized 
states and \1jj2) of the data qubit. 

(b) Elliptically polarized states: 

We also discriminate several pairs of elliptically polarized states. Experiments have been performed, using the pulse 
sequence given in Fig. (5), for a=90° and ellipticities e— 0°,15°, and 30°. For each value of e, we perform the 
experiment for 26»i=20°, 40°, 60°, 90°, 120°, 140°, 160°. As described above (in section III(i)), the data qubit states 
\ipi) or ^2) of Eqn. (5) are prepared respectively by applying a (2rj)(T r /2+4>) or a (2?7)_( 7r /2+0) pulse, where 77 and 4> 
are calculated from Eqn. (4). Ancilla qubit is prepared by using Eqn. (12). For a=90°, since \4> a)=\4)2) , ancilla qubit 
is prepared by (2?y)„( 7r /2+0) pulse. Figure (9) shows both experimental and simulated results of the probability of 
successful discrimination of pair of elliptically polarized states \ip\) and \ip2) as a function of 26\ (the angle between 
l^i) and |^2), as shown in fig. la), for different ellipticities for a fixed value of a=90°. From Fig. (9) one can obtain the 
probability of discrimination (P) of pair of elliptically polarized states, as a function of cllipticity. However for a=90° 
the maximum probability of discrimination for any ellipticity is always obtained for 26 \— 90°. The experimental 
results for low ellipticities (fig. 9) match well with the theoretical results, but deviates for higher ellipticities. Similar 
results have been obtained in o ptic s, where partial measurement in the Bell basis has been done for the discrimination 
of elliptically polarized states. [38j. 

IV. CONCLUSION 

The implementation of a programmable quantum state discriminator by NMR has been demonstrated. The device 
discriminates pair of data qubit states unambiguously (error free) that are symmetrically located around some fixed 
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state. One can use the same device (without changing it's parameters) to discriminate any pair of data qubit states, 
by suitably preparing the ancilla qubit. However the probability of discrimination depends on the parameter of the 
device (angle a). It may be noted that since NMR is an ensemble measurement, it is inevitable that to do projective 
measurement one has to prepare the input state twice. The probability of successful discrimination is obtained as a 
function of the angle between pair of data qubit states and the rotation angle of the unitary operator of the protocol. 
The states of the ancilla (programme) qubit that represent different programs can be nonorthogonal, which indicates 
the quantum nature of the programming. It is further shown that if the pair of data qubits are in elliptically polarized 
states then the probability of successful discrimination is also a function of ellipticity. 



Useful discussions with Arindam Ghosh and Karthick Kumar are gratefully acknowledged. The use of DRX-500 
NMR spectrometer funded by the Department of Science and Technology (DST), New Delhi, at the Sophisticated 
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Unitary operator corresponding to a radio frequency (r.f) pulse of angle a and phase (direction of r.f pulse) <f> is, 
Rj,(a), which is also called as (a)^ pulse, 



where n is a unit vector whose direction is along the direction of r.f pulse and / = I x i + I y j + I z k 1 where I is the 
angular momentum operator of spin 1/2 nuclei. 



and z-axis (direction of static magnetic field), and </> is the angle between h and x-axis. Here 9 — 90°, since r.f pulse 
is applied perpendicular to static magnetic field. 

After simplification, unitary operator of (a)^ pulse, R$(a) can be written as, 
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APPENDIX 



R<t,{a) = e 



—ian.I 



In spherical polar coordinates n=I x cos((p)sin(9)i + I y sin((j))sin(6)j + I z cos{9)k 1 where 9 is the angle between n 




-e"*sm(a/2) 
cos(a/2) 




10 



Here = 0° gives (a) x pulse, and 4> = 180° gives {a)- x pulse. Similarly <j> = 90° and = 270° gives (a) y and 
(a)- y pulses respectively. From R^(a), one can calculate any unitary operator, corresponding to any arbitrary angle 
and phase. For example the unitary operator corresponding to (277)^/2+0) pulse is, 

(cos(rj) — e~ l< ^ sin(rj) \ 
e lv sin{ri) cosyq) J 

The unitary operator of (27y)_( 7r / 2 +0) pulse is given by the Hermitian conjugate of the above. 
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FIGURE CAPTIONS : 

FIG. 1: (a) Pictorial representation of elliptically polarized states j^i) and \1jj2) of the data qubit. 
They are symmetrically placed with respect to |0). When 2# 1 =90° the two states are orthogonal. 
Ellipticity e is defined as, e = y/x. When y=0, e = 0°, and are linearly polarized states. 

(b) Pictorial representation of linearly polarized states j^i) and \ip 2 ) of the data qubit, \ipA) is the 
ancilla qubit. When the data qubits (l^i) and ^2)) are elliptically polarized states as shown in Fig. 
(la), then ancilla qubit \iPa) is also elliptically polarized state (not shown in fig. (la)). 

FIG. 2: Quantum circuit for the discrimination of data qubit state \4>d)= or |^ 2 ), using an 
ancilla qubit, prepared in state \iPa)- The unitary operator U needed for such a protocol consists of 
two CNOT gates, two NOT gates (X) and two other single qubit gates Ui and u 2 . For projective 
measurement controlled-<j 2 gate is needed at the end of the of the protocol, as described in the text. 

FIG. 3: The pulse sequence for creation of pseudopure state from the equilibrium state for a 
proton - carbon-13 two qubit system, using the method of spatial averaging |53[. In the product 
operator formalism |54|. equilibrium magnetization can be represented by AIi z +I 2z) where 1 stands 
for proton and 2 stands for carbon (since 71 # ~ 4713^). All the pulses are applied on 1 H so there 
is no change in carbon magnetization. 4/ l2 is converted to 2(/ l2 --y/(3)/i, / ) by (ir/3) x pulse, and 
gradient pulse kills the transverse magnetization 2^/(3)/i 3/ . The remaining magnetization of 1 H, 
2I lz is converted to y(2) ( I lz — Ix y ) , by a (tt/4:) x pulse. Evolution under J-coupling for time 1/2J (i.e. 
evolution under the unitary operator e~ l7Tllzl2z ) converts it to \J{2)(Ii z + 2Ii x I 2z ), which is converted 
to (Ii z — Ii x ) + (2I lx I 2z + 2Ii z I 2z ) by a (7r/4)_^ pulse. At the end a gradient pulse is applied to kill the 
transverse magnetization, yielding the magnetization (hz+hz+^hzhz) , which is a 1 00) pseudopure 
state. All the pulses are applied at resonance so that chemical shifts are refocused through out the 
pulse sequence. 

FIG. 4: (a) Equilibrium l H and 13 C spectra of 13 CHCl 3 dissolved in CDCl 3 . 

(b) spectra obtained after the preparation of pseudopure state by using the method of spatial 
averaging using the pulse sequence given in Fig. (3). To obtain these spectra, 7r/2 read pulses are 
used on each spin. The appearance of a single resonance line with positive intensity for each spin 
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(double the intensity of carbon and half that of proton compared to respective equilibrium spectra 
(fig. 4a)), is a confirmation of the 1 00) pseudo pure state (hz+hz+^hzhz)- 

FIG. 5: The pulse sequence for implementation of the quantum circuit of Fig. 2. The data 
qubit ( 1 if) is prepared in elliptically polarized states j^i) and \ip 2 ) (eqn. 5) by (2r))^/ 2+ ^ and 
(2?7)_( 7r /2+</>) pulses respectively and ancilla qubit ( 13 C) is prepared in state \ipA), by (2^)_( 7r / 2+( ^) 
pulse for a = 90°. In case of linearly polarized states (eqn. 2, y=0), data qubit is prepared in states 
l^i) and \tp 2 ) by (29i) y and (26 , 1 )_ J/ pulses respectively and ancilla qubit is prepared in state \ipA) by 
(26*2)2, pulse, where 202 is calculated according to Eqn. (12). Figure (2) contains four single qubit 
gates and two CNOT gates. NOT gate (represented by X in fig. 2, eqn. 9) is implemented by 
ti x pulse. Ui and u 2 (eqn. 9) are implemented by (a)- y and (a) y pulses respectively. CNOT gate 
(eqn. 10) is implemented by the pulse sequence /2)l-(ir /2) y -(l /2J)-(ir /2)1-(it /2) 2 _ z , where (vr/2)^ 
pulse is obtained by the composite pulse (n / 2) y -(7i / 2)l-(n / 2)\ y and (rr/2) 2 _ z pulse is obtained by the 
composite pulse (11 /2) 2 _ x -{iv /2) 2 -{ii /2) 2 X . The first (n/2) 2 _ x pulse of composite (-z) pulse is canceled 
with the last (n/2) x pulse of CNOT gate. All the pulses are applied at resonance, such that the 
chemical shifts are refocused throughout the pulse sequence. 

FIG. 6: Proton spectra of 13 CHCl 3 obtained after the implementation of pulse sequence given in 
Fig. (5), where the initial states of data and ancilla qubit are prepared in linearly polarized states 
(section III (i) ) , 

(a) |V>d)=|-0i>, 2^1=90°, and a=90° 

(b) |^£,)=|^i), 2^1=90°, and a=60° 

(c) |^£,)=|^i), 2^1=90°, and a=45° 

(d) |^£,)=|^i), 20!=9O°, and a=30° 

(e) \if> D )=\il> 2 ), 201=90°, and a=90° 

(f) |^ D )=|-0 2 ), 2^1=90°, and a=60° 

(g) \il> D )=\il> 2 ), 2^=90°, and a=45° 

(h) |^£>)=|^i), 2^1=90°, and a=30° 

In each of the above experiments \ipA) is initialized by choosing 29 2 to satisfy Eqn. (12). A 
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complete set of these experiments have been carried out for different values of a by varying 29 1 and 
26 2 (satisfying eqn. 12). The results are plotted in Fig.s (7,8). 

FIG. 7: Probability of discrimination P (resultant intensity of the transition of the data qubit, 1 H) 
from Fig. (6) and corresponding experiments for various 29 \ and a. The 29 2 is adjusted to satisfy 
Eqn. (12) in each case. The expected intensities (shown by thick lines) are obtained by simulation 
of the pulse programme of the pulse sequence given in Fig. (5) using MATLAB programme. Since 
the total pulse sequence lasts for about 11.8 ms, which is much less than 7\ and T 2 of both l H and 
13 C, the relaxation effects were not included in the simulation. However all the experimental data 
points are normalized with respect to the experimental spectrum of a = 90° and 29 1 = 90° for which 
the intensity is taken as 0.5, which is the theoretical expected intensity. The maximum probability 
of discrimination (P max ) is obtained for 26'i=90 o for all values of a. However, the value of P ma x 
depends on the value of a. 

FIG. 8: The results shown in Fig. (7), are replotted as a function of a for various 29\. The 
continuous curves are simulated results and the experimental data points are shown by crosses. 
From these curves one can find the optimum value of a for a given 29\. 

FIG. 9: Experimental and simulated results of probability of discrimination (P) of pair of 
elliptically polarized states j^i) and j^) are shown , as a function of ellipticity (e) and 29\ (angle 
between j^) and ^2)) f° r ct=90°. Simulated results sans relaxation are shown by thick lines. 
However all the experimental spectra are normalized to 29\ = 90° for e = 0° to the expected value 
of 0.5. 
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